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picture. 



PACS numbers: 03.65.-w, 02.30.Mv, 02.30.Tb, 02.60.-x, 02.60.Lj, 05.30.-d, 05.30.Jp, 
05.30.Fk 



Submitted to: J. Phys. A: Math. Gen. 



\ Postdoctoral fellow at the Doppler Institute for Mathematical Physics and Applied Mathematics, 
Prague, Czech Republic 



Boundary effect of a partition in a quantum well 



2 



1. Introduction 

Quantum singularity is a point defect in an otherwise 'regular' system treated in 
quantum mechanics. Despite the simple setting, such a system exhibits various features 
which are intriguing both mathematically and physically. On mathematical sides, it 
permits rich structures admitted under the general class of self-adjoint Hamiltonians 
describing the system. The mathematically allowed class is rather large, but when the 
system is linear (i.e., one dimensional), it is concisely characterized by the unitary group 
U(2) (see, e.g., [HOE]) which includes, among others, the familiar singular interactions 
realized by the Dirac delta potentials with arbitrary strengths. On physical sides, one 
the other hand, the quantum singularity admits interesting phenomena such as duality, 
super symmetry, anholonomy (Berry phase) and spontaneous symmetry breaking [U [5j [6] 
(see also [?]). An application for a qubit device, which might become realizable by future 
development of nanotechnology, has also been suggested [Sj. 

All the studies on quantum singularity mentioned above have assumed one-particle 
systems at zero temperature, but since actual realizations of the systems would involve 
many particles acting under finite temperatures, it should be inevitable to take into 
account the statistical aspects of the particles. The aim of this paper is to present a 
case study of such statistical aspects by analyzing a system which is both simple and 
familiar - a quantum well. Specifically, we consider particles in a quantum well with a 
partition placed in the centre. The partition is made by an infinitely thin wall which 
forms a single point defect inside the well. The partition separates the well into two half 
wells, each of which is supposed to contain iV identical particles. For the statistics of the 
particles, we consider separately the two cases, the Bose-Einstein and the Fermi-Dirac 
statistics, and the entire well is put under various temperatures, ranging from zero to 
infinity. 

To characterize the partition in a simplest nontrivial term, we assume that the 
partition enforces the Dirichlet condition on the left and the Neumann on the right 
(see figured]), which is perhaps a simplest nontrivial combination permitted quantum 
mechanically. We focus on the pressure, or statistical quantum force, acting on the 
partition, which arises from the discordance in energy levels between the two half wells 
[9]. We shall be content to work in the one dimensional setup, but the extension to multi- 
dimensions will be done analogously by considering the class of singularities pertinent 
to the dimensions. 

The temperature behaviour of the net quantum force — the difference between the 
pressures from the left and the right of the partition — exhibits a number of interesting 
and unexpected features. For instance, the force curve as a function of temperature 
starts off with a finite value at the zero temperature and reaches a single minimal point 
at a finite temperature before diverging in the high temperature limit according to a 
certain scaling law, i.e., the square root of the temperature. This overall temperature 
behaviour can be observed for both the bosonic and the fermionic cases, but in the latter 
case the curve in the low temperature regime shows an additional little twist which is 
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Figure 1. Quantum well with a partition at the centre. Under the Dirichlet (ip = 0) 
and Neumann (ip' = 0) boundary conditions on the left and the right of the partition 
(while we choose the Dirichlet for the two ends of the well) , the eigenstates (the lowest 
two are shown in both half wells W ) have different energy levels. When the same 
number N of particles are introduced in each of the half wells, these level differences 
lead to a net force on the partition which exhibits interesting temperature behaviours. 



absent in the former. Another feature worth mentioning is the distinct characteristics 
of the dependence on the number of particles; e.g., both the zero temperature value and 
the minimal force are of the order of N for bosons while they are of the order of N 2 for 
fermions. 

Our analysis on the quantum force is carried out both numerically and analytically. 
Some of the initial results have been reported earlier [TOl ITT] , and here we shall 
present the full detail of our complete analysis, including the analytical approximations 
performed with assistance of numerical solutions of some transcendental equation. We 
shall see that our goal to find analytic formulae that can account for the pure numerical 
results as well as the statistical feature of the particle number dependence mentioned 
above is achieved reasonably well. Our outcomes illustrate how quantum singularity, 
realized by a partition in the present model of a quantum well, can give rise to physically 
measurable effects in a plain form such as pressure. 

This paper is organized as follows. In section [2] we provide the basic account of 
our model as well as the overall features of the quantum force obtained by numerical 
computations for all temperature regions. In section [31 we present the analysis for the 
high temperature limit where the force exhibits the common scaling behaviour for the 
two types of particle statistics. Section H] is devoted to the analysis of the medium 
temperature regime where the force curve has a minimal point. Our argument for the 
analytic approximation is given for bosons and fermions, separately, since the statistical 
difference becomes important in this regime. Similarly, section [5] is devoted to the 
analysis of the low temperature regime where the force possesses the characteristics of 
the iV-dependence and further develops the twist in the fermionic case. In section [6] 
we briefly mention two physical quantities, the shift of partition and the transfer 
of particles between the two half wells, as alternatives which may be more directly 
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observed than the quantum force. Finally, we present our conclusion and discussions in 



section [7J Appendix Appendix A gives the outline of our numerical computation, and 
Appendix Appendix B contains the technical detail of our approximation employed in 
section @J 

2. Quantum force on a partition and its temperature behaviours 

Before we present our full analysis of the quantum force, we provide the basic account 
of the system for which we discuss the boundary effect for various temperature regions. 

The system we consider is a one- dimensional quantum well W possessing a partition 
wall at the centre, given by the interval [— I, I] with the partition at x = 0. To define the 
system in quantum mechanics, we need to specify the boundary conditions for the wave 
function if)(x) imposed both at the two ends of the well x = ±Z and at the partition on 
its left and right sides x = ±0. In quantum mechanics, allowed boundary conditions 
are those which respect the probability conservation requirement. If we assume that 
the partition is impenetrable, i.e., the particles cannot penetrate the partition, then the 
most general form of the boundary conditions is given by 

V(±Z) + L ±l tfj'(±l) = 0, ^(±0) + L ±0 V'(±0) = 0, (1) 

where L±i, L±q are arbitrary real constants including infinity. For instance, L + q = in 
(pQ) implies the Dirichlet condition ip(+0) = whereas L +0 = oo implies the Neumann 
condition ip'(+0) = at x = +0. Similar choices of boundary conditions should be 
made at the other three points, x — ±Z and x = — by specifying the constants L±i and 

In this paper, to realize a nontrivial partition in the simplest setting, we adopt the 
set of boundary conditions provided by L± t = 0, L +0 = oo and L_ = 0, that is, the 
Neumann boundary condition at the right side x = +0 and the Dirichlet at the left side 
x = — of the partition as well as at the both ends x — ±Z of the wall: 

VH) = o, <M-o) = o, <//(+o) = o, v(0 = o. (2) 

The quantum well W is then split into two half, physically distinct, wells W + and W~ 
by the partition. In the two half wells W with the boundary conditions ([2]), the free 
Hamiltonian 

h 2 d 2 

admits energy levels 



H = ————— (3) 

2m dx 2 



h 2 /7T\ 2 

E n = e n £, £ = — (-), n = 1,2,3..., (4) 



2m V I 

§ The allowed set of boundary conditions |H) forms the group [C/(l)] 4 representing the possible 
self-adjoint domains of the free Hamiltonian in the interval [— I, I] divided into two. If we allow the 
partition to transmit particles, then the boundary condition can be more general than ([T]) and given 
by U(2) x [C/(l)] 2 (see Q] GOT- 
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Figure 2. Illustration of the distribution of particles over the levels in the two half 
wells at a low temperature. In the fermionic example s = 1/2 is assumed. The gap of 
the levels and the difference in the populations in the two wells induce a non- vanishing 
force on the partition, which is dependent on the temperature as well as the particle 
statistics. 



which are distinct for the two half wells (see figured]), 

(n-\) 2 iorW\ 
n 2 for W . 

It is then expected that the level gap between the two half wells W ± gives rise to 
noticeable physical effects, and one obvious example will be the force (or pressure) 
acting on the partition due to the difference of the pressures exerted by particles inside 
the two half wells. To extract the pure boundary effect, we put an identical number of 
particles in each of the half wells, and thereby study the temperature dependence of the 
net force emerging on the partition. We do this for the two cases of particle statistics, 
bosons and fermions, which possess different features because of the different statistical 
distributions over the energy levels (see figure [2]) . 

To be more explicit, let N denote the number of the identical particles in both of 
the two half wells W in each of the 2s + 1 spin degrees of freedom, where s is the spin 
of the particles. At temperature T, the particles distribute over the levels according to 
the population number (per spin degree of freedom) 

N n = TT~~ , (6) 

n e a+be n _ ' v I 

where we have introduced the statistical index, 

,_ 1)2s f 1 for bosons 
| —1 for fermions 

and the shorthand 

8 



with ks being the Boltzmann constant and 8 given in (j3J). The temperature-dependent 
quantity a is determined by the total particle number constraint, 

iV = £iV n . (9) 
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The force acting on the partition from one side is then given by 

F = -(2s + 1) E N n — n - = (2s + 1)— E N n e n , (10) 

n L n 

which is simplified as 

/ = £iV»e„, (11) 

n 

in terms of the reduced (dimensionless) force defined by 

/ = — ^F. (12) 

J 2s + 12S K J 

Denoting by f ± the forces on the partition in the half wells W ± , we have the net force 

on the partition (from the left to the right), 

A/ = r - / + = £ Ken ~ £ N >t- (13) 

n n 

Throughout the paper we use the superscripts ± to refer to quantities pertaining to the 
half wells W ± , respectively, and also use A to refer to the difference between the two 
quantities such as (|T3|) . However, we will often omit these superscripts for brevity as 
long as the distinction is unnecessary in the argument. 

Our aim will be to determine the net force Af in ( TT3l) as a function of the 
temperature. For our later convenience, we introduce the reduced temperature t by 

( = I = |t. (14, 

Numerical evaluations of Af in terms of t can be achieved readily for finite N, and 
the outcome for N = 100 is shown in figure [3] for the whole range of t. We observe 
immediately there that (i) the net force has a finite zero temperature limit, (ii) it 
exhibits a power law (as it approaches the dotted straight line) in the high temperature 
limit, and (iii) it has a single minimum in between. These overall features are common 
both in the bosonic and fermionic cases, and can be seen for other values of N different 
from N = 100 as well. The basic difference between the bosonic and fermionic cases 
lies in the fact that the strength of the force is of the order of N for bosons while it is 
of the order of iV 2 for fermions. Interestingly, the same difference in the N dependence 
can also be seen for the temperature value of the minimum force. In most part of our 
discussions, the numerical results will be presented for N = 100 particles in each spin 
degree of freedom (which is already a realistic population number in nanoscale quantum 
experiments [12]), but the conclusion remains the same for larger N. 

Analytical evaluations of the force, in contrast, are not readily obtained, since 
the required sums in (fT2j) cannot be performed exactly, and an analytical solution for 
a which has to be determined is also difficult to find. We are thus forced to apply 
some approximations, but the standard approach of approximating the sum with an 
integral turns out to be insufficient to get the accuracy needed in the present situation, 
leaving us to work out novel methods for the present problem. For instance, the Fermi- 
Dirac integral, which is often used to approximate the sum and is related to the Lerch 




Figure 3. Double logarithmic plot of the net force Af as the function of the 
temperature variable t, for bosons (left) and fermions (right), at N = 100, obtained 
by a numerical computation (solid line) . Their common high-temperature asymptotics 
([23]) is also displayed (dashed line). 



transcendent, can be expressed via an asymptotic series [13], but the truncations of 
the series do not provide enough precision to recover the force difference obtained 
numerically. As we shall see soon, our approximations developed in this paper and 
performed partly with the help of numerical solutions for transcendental equations 
provide formulae which fit the numerical results reasonably well. In particular, the 
three salient features (i), (ii) and (iii) mentioned above will be seen to be all reproduced 
properly by our analytic methods. 



3. High temperature regime 

As we have seen from the numerical results mentioned in figure [3j in the high- 
temperature limit the force Af exhibits a certain scaling behaviour which is common to 
both the bosonic and fermionic cases. We first show that this scaling behaviour can be 
explained analytically based on a rather simple argument which is completely analogous 
in the two cases. 

To study the high-temperature regime, we first note that the population number 
Ni decreases for increasing temperature (see (jBJ) for n = 1) and, accordingly, we expect 
a to increase to higher positive values for large t (or small b). It follows that the factor 

q:=e~ a (15) 

will be extremely small in the high-temperature regime (recall that a is a function of 
t), and this leads us to expand N n in q as 
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with rj given in (171), which is valid for any positive a. Thus we find 



-kbe n 



=E(w)' 
fe=i 



a 
2 



1 

2 



E e 

!= — OO 



-kbe r 



(17) 



»^ = i7Etf»=E(w)*E e " 

n=l fc=l n=l 

where we have introduced the constants a + = , o~~ = 1 corresponding to the half 
wells W , and extended the meaning of the notation e n [see (JSJ)] to negative n as well. 
Applying the Poisson summation formula 



OO OO ™ 

e = e y d ^i/( w ) e 



2irimu 



we obtain 



^ = E(w)* 



fe=i 



2 



4t E (2^-1' 



e kb 



(19) 



where we note that the factor 2<7 ± 
find 



1 is simply =pl. Similarly, for the force /, one can 



/ = E N ^n 

n=l 

oo 



k=l 



71 



lQk 3 b 3 



E (2a 1 



2^ 



-m 



2 „-?rT 



(20) 



In the high-temperature asymptotic limit 6 — > 0, we have q — > and hence it 
suffices to consider only the first few terms in the sums over k [both in ( fl9|) and ( |20l) ]. 
and within each term to keep only the m = term in the sums over m (the m/fl terms 
being exponentially suppressed). Now, the leading k = 1 term in (fl9l) gives 



1/2 



2iV 



7T 



+ 0(6) 



(21) 



This result shows that, for high temperatures, a tends to infinity logarithmically. Since 
this leading behaviour of q is independent of a, inserting it into ( 1201) gives that the 
leading 0(b~ l ) term of / (coming from k — 1, m = 0) is also cr-independent. Hence, 
the contribution coming from this term will cancel out between the two forces / + and 
/~ in the net force f|T3|) . 

The first nonvanishing contribution for the net force comes from the first subleading 
term in q. Incorporating the k = 2 term as well for q, we find 



2N 



b_ 

71 



1/2 



+ 2N |er — r]V2N] - + (& 3/2 ) . 



Plugging this into (120]) and calculating the net force f|T3l) we obtain 



A/ 



N 

T 



7T 



+ it 



(22) 



(23) 



which shows that the net force diverges asymptotically as t —>■ oo according to the 
square root of t, and that it is proportional to the particle number N. Note that these 
are true both for the bosonic case and the fermionic case. These results explain the high 
temperature asymptotic behaviours observed by the numerical analysis in figur^Sl 
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By incorporating higher orders, the approximation may be improved easily. For 
example, by taking account of the next order, one finds 



N ft\ 1 ' 2 



7T 



N 

71 



(V2- 



i) V N - i 



(24) 



This improvement, however, is not sufficient for describing the behaviour of the quantum 
force in lower (medium) temperature regimes where the force takes its minimum. 



4. Medium temperature regime 

Apart from the strength of the force, the overall shape of the force curve - for all 
temperature regimes - agrees for bosons and fermions. Indeed, for both of the statistical 
cases, the force starts with a finite value at t — followed by a plateau and then 
by a steady decrease for larger t, and it exhibits even quantitatively the same high- 
temperature behaviour for t — > oo as discussed in section 3. Moreover, the force admits 
only one minimum between the two limiting domains (see figure [3]) for the two cases. A 
closer inspection of the curve reveals, however, that in the medium temperature regime 
the two cases differ in scalings with respect to the particle number N. Explicitly, our 
numerical results (see figure H]) show that, for AT 3> 1 , the value of the minimum A/ m j n 
and the temperature £ m i n where the minimum occurs are both proportional to N for 
bosons, while they are proportional to A^ 2 for fermions. We now give an account of these 
distinctive features characterizing the statistics of the particles by analytic means. 



4-1. Bosons 

To analyze the regime of temperatures including the minimum t min for the bosonic 
case, we note that interpolating between the slow logarithmic increase of a for high 
temperature and the low-temperature behaviour fl64|) of a (which will be derived in 
section [5]) suggests the existence of a temperature regime fulfilling 

6<1, |a|«l. (25) 

This is the 'medium temperature regime' we wish to consider here. Our procedure to 
derive the force difference is based on an approximation for a fulfilling these conditions 
in (l2"5l) . and we confirm the consistency of our argument by examining the validity of 
the conditions later. 

To proceed, we first solve the constraint ([9]) for a by approximating the sum as 

N = E —ur — t « E — -r- ■ ( 26 ) 



n=l 



C 



a +be n _i + be r 



This approximation must be good at least for the lower levels — the ones that provide 
the dominant contribution in the sum (note that N n falls rapidly with n) . Let us write 
this in the rearranged form 

N ™ 1 

— « E ■ ( 27 ) 
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Figure 4. (Left) The minimum of the force curve, for bosons, at N = 100. Solid 
line: numerical computation; dashed line: the approximation (|40[) : dotted line: the 
approximation (|44p . (Right) Numerically observed scaling behaviours of the minimum 
force A/ m ; n and the temperature £ m i n at which the minimum occurs. The limiting 
values determined in ([55)) are also displayed. 



If one solves this in favor of ta, then the solution will be a function of 4f. Using the 
standard expansion formula (see, e.g., page 42 of [H]), we find that the summation in 
( 127)) is evaluated as 

oo i r 7rtanh7rv / ta £ Qr 

S(ta):=V - = ^ ' (28) 

where we have used 

{ - for 

r := I 2 ' (29) 

1 for IV". V ; 



Both formulas in (128)) are valid for negative a as well, and are a smooth function of ta at 
zero, see figure El Then, for a given t, or given t/iV, we should solve the transcendental 
equations 

S ± (ta ± ) = — (30) 

to obtain ta + and ta~ which are the values corresponding to the half wells W + and 
W~ , respectively. 

At this stage, it is informative to determine those temperatures at which either a + 
or a~ vanishes. Let t§ and t§ be the temperatures such that a + (to) = a ~(to) — 
hold. Using ( 1281) . one can immediately determine t$ as 
+ _2N __QN 

t - t - —. (31) 

One then finds, numerically, 

t+a~(t+) = -0.7627, a + (t ) = 0.9026 . (32) 
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Figure 5. The functions S + (dotted line) and S~ (solid line) as defined in (|28|) . 
They are smooth at zero, and diverge at —ef resp. — where the n — 1 term in (|28p 
diverges. The dashed line depicts the approximation for S + around 1 that is based on 

(El. 



Our assumptions b <C 1 and {a^l -C 1 are in fact valid at these temperatures, since we 
have ~ N ^> 1 and olq ~ l/t — b ~ 1/iV <C 1. Obviously, the same will apply for 
the temperatures between and in the neighborhood of these two values . Moreover, 
we can expect from these results that, for t ~ N in general, 

o:± ~ b ~ 1/JV< 1. (33) 

In what follows, we shall be interested in this temperature regime, where the minimum 
of the force difference A/ min is found in the numerical result (see figure H]). One can 
show — see Appendix Appendix B for the details — that, for t ~ N , the relative error 
of our approximation ( 1261) decreases with increasing N. 

Next, we derive an approximated analytical expression for the net force Af obtained 
under (1331) . We do this with the help of an integral approximation for the force sums. 



Presenting all the technical details of the calculation in Appendix Appendix B here we 
summarize only the result, 

-A(ta). (34) 

This formula has a relative error that tends to vanish in the large- N limit. 

We observe from (134j) that, in this medium temperature domain, Af /N is also a 
function of t/N as tat is. Consequently, the force curve Af(t) is scale invariant in iV 
(i.e., it preserves its shape under the rescaling of N) as both the effective range of t 
and Af scale linearly with N. As a special case, we see that the temperature of the 
minimum force occurs at t m j n ~ N, A/ m j n ~ N, confirming our observation made in the 
numerical analysis in figure HI 

Let us now determine the ratios t/N and Af /N in the large- N limit. One possible 
way to do this is to use various values of t/N to solve numerically ( |30l) for a and insert 
the outcomes into (1341) . and thereby reproduce the force curve at those values. For 
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example, the minimum of Af /N as the function of t/N can be estimated in this way, 
yielding the ratios 

^^ 0.5480, 0.5967. (35) 

These values agree with the numerical results pretty well as readily confirmed in figure HI 
To find an analytical formula for Af as well, one may use some expansion 
approximation of the functions S ± [see (1281) ] around some chosen value. As an example, 
let us choose t to be in the vicinity of . There, \ta~\ <C 1 , and the equation we wish 
to solve becomes 

N irV tor coth nV tor — 1 
T ~~ 2ta- 

6 90 945 V ; ^ V ' > K ' 

Inverting this, one can obtain the solution as 

5/t 6\ 5n 2 / t 6\ 2 , n (ft 6\ 3 \ 

As for to + , we observe that it is close to the value 0.9 [see (132]) ]. where irVta + ~ 3 and 
tanh7rv / to+ already almost saturates to its large- variable asymptotic value, 1. Taking 
this simple asymptotic approximation, 

tanh7r\/to+ » 1, (38) 



one can determine tot + by the solution of N/t = 7r/(2vto+) , that is, 

7T 2 / t ^ 2 



-l--- 2 ) +-S- (4°) 



Applying the solutions (1371) and (139]) in (|34|) yields 

Af ^ 7T 2 / t 6 \ 2 6 
~~N ~ 14 ViV ~ ^ J + n 
This quadratic formula (HUT) indicates that the location of the minimum occurs at 

^ = 0.6079.... (41) 
N N n 2 K ' 

Figure H] shows that around the minimum the formula (I40p reproduces the numerical 

curve of the force reasonably well. 

To improve our formula to achieve a better agreement between PI"]) and (133]) . we 

consider the approximation: 

tanhx* + (x - x*) 

tanh x w ^ ^- . (42) 

1 + tanhx* • [x — x*) 

This formula is precise up to the quadratic Taylor term at x*, and behaving much better 
than the quadratic Taylor polynomial approximation in a larger neighborhood. For the 



expansion point x* = it J (ta + )* , we may simply choose the above-mentioned value, 3. 
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Using the notation t* /N := S + ((ta + )*) — which is close to t /N — and assuming an 
expansion of the form 

to+ - < to+ >* = * {jf-v) + « (jt ~ ^ + ° ((* " v)) • < 43 > 

we obtain a quadratic approximation for ta + that is better than (139]) . Accordingly, the 
force difference is obtained in the improved form, 
Af ft \ 2 

— f « 0.5121 • — - 0.5465 +0.5967, (44) 



N \N 
having its minimum at 

^^ 0.5465, ^Jp« .5967. (45) 

We remark that the parabolic approximation used above is valid only locally and not 
suitable for describing the whole medium temperature domain. To find a formula valid 
for other values of temperatures, we may simply expand at the value of the interest 
and/or using approximate formulas for the other parts of the functions S. A more 
universally valid approximate formula, which describes the force curve in the whole 
medium temperature region, is obviously desirable but it is rather difficult to find at 
the moment. 



4-2. Fermions 

The numerical results shown in figure H] suggest that, for fermions, the medium 
temperature regime that includes the temperature of the minimum force may be defined 
as the domain where t ~ N 2 , Af ~ N 2 are fulfilled. Now we investigate this regime to 
seek an analytic approximation of the force curve there. 

As we did for the bosonic case, we first determine a. Our integral approximation 
presented in Appendix Appendix B gives 

w = |5=STT M ^T + ^ I(o) (46) 

with the Fermi-Dirac integral 1(a) given by 

I( a ) ;= r d f . (47) 
v ' Jo e a+ y +1 v ' 

Since e <t +l ~ N\ < 1 , for large N we can simplify ( f46l) to 

Ija) 

where, for temporary convenience, we solve our condition for i as a function of a, 
instead of the reverse. We observe that t ~ iV 2 implies a = 0(N°), unlike in the 
medium temperature regime of the bosonic case. 

Since there arises no difference between a + and a~ in this leading order, we need 
to consider the subleading term for Aa. This can be derived from fj4*6|) as 



1- iV 

N m Vt 1(a) or t w — -rr , (4* 



= AiV ps A 



e a + l 



4=/ / (a)Aa, Aa^— - , (49) 

Vb V ; 2N(e a + l)I'(a) V ; 
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where we have used the fact that, in leading order, the difference between a + or a~ can 
be ignored (and we have again eliminated i as a function of a) . 
For the force, our integral approximation provides 

f = § 1 = 6 ~ 3/2 fj e^vl + 1 ~ °^)+ el2 j e *ly* V +1 i ( 5 °) 

and hence by integration by parts and (|49p we obtain 



with 



A/ w t 3/2 / s y 2 dy = Aa / — w — J (a , 51 



o 



— 1 —2 
J(a) := (e° + l)/(a)/'(a) = (e a + 1) d/ 2 (a)/da ' (52) 
All these approximations will improve for larger values of t and iV. We can see in 
(IBTj) that, having a = O(iV ), not only t but also Af scale with iV 2 in the medium 
temperature regime, and consequently the force curve is again scale invariant in N 
(shape-preserving under rescaling N) in this regime. 

As before, for a given t, one can obtain a and Af by solving a transcendental 
equation numerically or by using some analytical approximate formula. For instance, 
determining the location of the minimum — in essence, the minimum of J(a) — via 
numerical solution, one finds 

a min = -2.567, J(a min ) w 1.813, t min w 0.4447V 2 , A/ min « 0.453iV 2 . (53) 

These values are in apparent accord with the fully numerical results in figure HI 

For an analytical approach, first one can take the standard asymptotic series for 
the Fermi-Dirac integral 1(a). According to [13J, in the interval a G [—3.696, —1.314] 
the truncation 

7T 2 1 



1(a) v— ol 



24 (-a) 2 



(54) 



of the asymptotic series is the best available approximation for 1(a). Unfortunately, 
this is not sufficient for our purpose because the approximated J(a) obtained from this 
does not possess a minimum. However, noticing that, for I' (a), a better approximation 
is obtained by the further truncated 

I'(a) « —L= (55) 
2v— a 



we can find, using this in J (a), a minimum in the temperature regime we are considering. 
By numerically solving the arising transcendental equation, one finds 

a min ^-1.95, J(cw) « 1.96, t min ^ 0.61LA 2 , Af min ^ 0.49iV 2 , (56) 

which are not quite precise compared to f l53lh 
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Figure 6. (Left) The minimum of the force curve for fermions, at N = 100 . Solid 
line: numerical computation; dotted line: the approximation corresponding to (|56[) : 
dashed line: the approximation corresponding to (|60[) - (Right) the scaling behaviour 
of t m - m and A/ m j n . The constants determined at ([53)) are also displayed. 



To find a second, better approximation of 1(a), we utilize the fact that its integrand, 
e a+ v $ +1 * s c l° se to a tangent hyperbolic function (reflected and shifted) on y G [0, oo) , 
for negative a. Specifically, we may write 

f°° P r / x-i 

1(a) ~ / - 1 — tanh v — a ( y — v — a ) dy 
Jo 2 L v n 

= — ^= \-2a + ln(l + e 2a )} (57) 
2y/—a 1 J 

with p := ■ Here, the integrand is chosen to be simple enough but still to reproduce 

the true integrand exactly at y = and decreases to the half of the y = value around 
the same point with similar steepness. Omitting 0(e 4a ) terms can reduce the obtained 
formula to 

1(a) w — (-2a + e 2a ) , I 2 (a) ^ p 2 (-a + e 2a ) . (58) 
2\J —a v ' v 7 

Since no further available simplification can render the resulting transcendental equation 
analytically solvable for a, we shall use f )58|) to expand our J(a) approximation around, 
say, a* = —2.5 to second order, and determine the minimum of the quadratic Taylor 
polynomial, which polynomial we can rewrite in the form 

J (a) ss 0.134(a + 2.48) 2 + 1.64. (59) 

This way we reach 

a min S3 -2.48, J(a min ) w 1.64, t min « 0.466iV 2 , A/ min ss 0.411iV 2 , (60) 

which agree better with the numerical results; see (153]) and figure HI Should one need a 
further improved formula, the approach can be repeated with some enhanced or more 
diligent choice for approximating the integrand in 1(a). 
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5. Low temperature regime 

Having analyzed the quantum force on the partition in high and medium temperature 
regimes, we now consider the force under low temperatures in this section. 



5.1. Bosons 

To discuss the bosonic case, let us first approach from the zero-temperature end. Since 
at exactly zero temperature all the particles sit on the lowest available level which is 
the ground state for bosons, we immediately obtain 

N 

/+(0) = efN = — , f-(Q)=e- 1 N = N. (61) 

This implies 

A/(0) = ^V, (62) 

which is nonzero and is proportional to N. 

When the temperature is slightly above zero, the particles occupy the lower excited 
levels in addition to the ground state, and the transition from the ground level to the 
upper levels is more extensive in the half well W + than W~ , because the subsequent 
energy levels have a smaller difference in W + than W~ (see (jSJ)). Consequently, the force 
difference will decrease as t grows from zero. Indeed, in the two-level approximation for 
bosons, where the higher levels are treated as still completely unoccupied, the net force 
is found to be |10| 



Af{t) « +(3e~ 3/t -2e~ 2/t ), (63) 

which accounts for the decrease, see figure [7J Note that Af(t) starts to decrease for 
t ~ 1 , irrespective of the particle number N. Incidentally, we mention that the low 
temperature behaviour for a is 

a » -be x + In ( 1 + -M , (64) 



N 

which is a straightforward consequence of the approximation N\ « iV 



5.2. Fermions 

To study the zero temperature limit for fermions, we recall that the lowest N levels are 
occupied at t — 0, and from this we obtain 

m ^ t ;^, r( o) = fx = g±ig±i> , (65) 

n=l n=l 

which implies 

Observe that, in contrast to the bosonic case, the force difference in the fermionic case 
is (for N ^> 1) proportional to N 2 in the limit t — > 0. 
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Figure 7. The low-temperature behaviour of the net force A/, for N = 100 . 
(Left) the bosonic case. Solid line: numerical computation, dashed line: the two-level 
approximation (|63|) . (Right) the numerical results for the fermionic case. Observe the 
stcplike shape, characterized by the temperature values at the two points of inflection. 
The inset displays t-^gg/N and t^^/N , as the function of N. They appear to tend to 
constants that are approximately 0.237 and 0.472, respectively. Within an error, the 
second number is just the double of the first one. 




We have learned that the fermionic net force A/(i) differs from the bosonic one 
quantitatively in the zero temperature limit A/(0) - - it is proportional to N 2 for 
fermions while the order is iV for bosons. Another quantitative difference can be 
observed in the first 'turning point', namely the temperature where A/(t) starts to 
decrease when the temperature t is increased from zero - the turning point occurs at 
around t ~ N for fermions (as confirmed numerically) while it occurs at around t ~ 1 for 
bosons. Besides, there is a qualitative difference between the fermionic and the bosonic 
cases - the fermionic curve exhibits a single, very small but unmistakable 'steplike' 
pattern (or a depression) for any N during the initial decrease as shown in figure [SJ 
This step occurs at a temperature proportional to iV (see figure [8]), and the net force 
differs there from the O {N 2 ) zero temperature value by only an O (N°) = O (1) amount. 

These numerically observed properties can be understood analytically as follows. 
For fermions, the analogue of the bosonic two-level approximation corresponds to the 
situation where only the occupation of the iVth and (iV+ l)th levels {i.e., the two levels 
closest to the Fermi level) are different from the zero temperature value. On each side, 
this imposes the approximate equation 

N N + N N+1 = 1, -Wjv+i — — it~~ - = 1-N N = 1 -,(67) 

^ ^ e a+be N+1 _|_ I e -a-be N _|_ \ V > 

which can be solved for a as 

a ~ ~\ ( e N + e N +i) ■ (68) 
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Let us observe that, with this approximation for a, 

NN + l = e «+bel +l + I ~ i-^W-l = e -a- feej ! +1 _, + J ~ e (2H)L + 1 ( 69 ) 

(/ = 1,2,...) in the leading order of N in the exponents. Keeping only two nontrivial 
levels for calculating the forces as well (i.e., I = 1 only), we have 

00 N N oo 

/= £eJV n = ]Te n -]Te n (l-Ay + ]T e n N n 

n=l n=l n=l n=N+l 

*f(0) + (e N+1 -e N )N N+1 . (70) 

At this point, we observe that a cancellation in the leading order of N takes place, 
since we have ~ ew+i ~ N 2 but ejv+i — e^ = 2N + (1 — 2r) (and, in general, 
ejv+; - ejv+i-i = (2/ - l)[2iV + (1 - 2r)] ) with r defined in flU}, which is only O(N). 
An even higher cancellation will occur in the force difference, as this leading O (N) 
difference is the same on the two sides and, according to (1691) . N N+ i is also the same on 
the two sides in the leading order of N. To obtain a nonvanishing contribution, we use 
for iV^ +1 the first-order Taylor approximation 

1 1 e x 



e x+Ax + 1 e x + 1 ( e x + x f 
with X = bN and appropriate AX ± . This gives 



AX (71) 



1 / 1\ e bN 



e bN + 1 V 2/ (e fe ^ + l) 2 

1 e bN 

w A/(0) + -rr; ? . (72) 

M ; e biV + l (e^ + 1) 2 v ; 

Now we can see that this last expression contains t only in the combination bN = N/t, 
which explains why the decrease of the net force from the zero temperature value starts 
at t ~ N . Further, it is also visible that this decrease at t ~ iV is only a small 
O(N ) =0(1) phenomenon. 

This two-level approximation is not sufficient to explain the step in the curve 
pertinent to the fermionic case (see figure ED- However, we may apply a semi- four- 
level approximation, that is, we assume four nontrivially occupied levels but use the 
two- level-approximated a. The net force evaluated by an analogous procedure then 
becomes 

1 3 bN e bN 13bN e 3bN 

~ A/ (0) + ^TT + ^ + 1 " ( e *w + l) 2 " (e«^ + l) 2 ' ( } 
From this, we are able to read off tbeg and t en d that characterize the location of the step 
in figure [8] by numerically determining the points of inflection as a function of t/N. The 
results 

t heg /N w 0.239 , W/iV w 0.426 (74) 
are indeed close to the numerically observable values 0.237 and 0.472, respectively. 
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Figure 8. The two- and semi-four-level approximations (dotted and dashed lines, 
respectively) near the fermionic step (solid line) for N = 100 . 



However, as seen in figure [HJ for the curve over the entire low temperature regime, 
this semi-four-level approximation is less satisfactory. Our further investigation shows 
that the net force Af is extremely sensitive to the error in a, and that incorporating 
two more levels provides only a smaller contribution compared to the change caused by 
the error. It seems, therefore, that a full four-level approximation for a is required for 
the net force to reproduce the numerical curve more precisely. 

6. The partition shift Al and the transfer number AN 

We have seen that the quantum effect caused by a set of nontrivial boundary conditions 
manifests in the force that arises on the partition in a potential well. For actual 
observation of the effect, however, there may be other quantities which are more readily 
measurable than the force itself. In this section we mention briefly two examples of 
such quantities: the shift Al in the position of the partition and the transfer AiV of 
particles between the two half wells (see figured]). For simplicity, our discussions are 
mostly restricted to the zero temperature limit t — > 0. 

To discuss the first example, suppose that the partition is allowed to move freely in 
the well. Due to the pressure, the partition will then move and acoordingly the widths 
of the half wells W change as 



so that the net force vanishes, AF(0) = 0. The portion £ of the shift Al = is thus 
determined by 



W + : Z — > Z(l — 



W~ : Z-Z(l+0, 



(75) 




(76) 
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Figure 9. (Left) the shift Al = It; in the position of the partition. (Right) the 
redistribution of particles AiV = N + — N. Both Al and AiV are determined from the 
stability, AF(0) = 0. 



From this we find 

1/3 

(77) 



j. r — 1 

f = with 

r + 1 



r(o) 



./ + (°). 

For the bosonic case, (162|) implies r = \/A and hence 

xfl- 1 

i = 4^ « 0.227. (78) 

This shows that the shift is rather large reaching nearly a quarter of the original width. 
Notice that for the bosonic case the portion £ of the shift is independent of the particle 
number N. 

For the fermionic case, in contrast, the force limit (1661) implies the factor r w 
1 + 1/(2 AT) for iV ^> 1 to each up and down spin, and hence the portion £ in (1771) 
becomes 

« K iF (79) 

which is quite small for large iV (e.g., less than 1 percent even for iV = 100). The result 
(179]) indicates that the fermionic shift is much smaller than the bosonic one (1781) and 
is almost undetectable for large N. We thus learn that the spread of particles over the 
levels according to the Fermi-Dirac statistics has the effect of balancing the partition 
near the centre. 

When the temperature t increases from zero, the shift £ decreases steadily for all t 
in both bosonic and fermionic cases. This is due to the fact that as t becomes higher 
the forces /~ and / + grow faster than their difference, rendering the ratio r closer to 
unity. This is seen clearly in our numerical analysis shown in figure [10J Interestingly, 
it also shows that, for fermions for instance, the turning point of the shift £ after the 
initial plateau in the vicinity of t = has a curious scaling property under the change of 
the particle number N. Namely, the inset of figure [10] indicates that the turning point 
temperature of £ depends linearly on iV (since, for iV£, it scales as iV 2 ), which is the 
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Figure 10. The cquilibrial shift of the partition as the funtion of temperature 
obtained numerically for N — 100; (Left) for bosons and (Right) for fcrmions. The inset 
in the fermionic case shows the quadratic iV-dependence of the combination N^it). 



same scaling law found for the force Af for its turning point in the low temperature 
regime. 

In our second example, we suppose that, unlike in the previous situation, the 
partition stays at the centre but instead we consider to move particles from one half 
well to the other until the equilibrium AF(0) = is achieved. Let N + and AL be the 
redistributed number of particles in the half wells W ± , i.e., 

W + : iV — > N + , W~ : iV — > AL, with N+ + N- = 2N. (80) 

The amount AiV = N + — N of moved particles with respect to the original numbers iV 
is another good measure of how far the original situation is away from the mechanical 
equilibrium. The condition to determine this number AiV is just r = 1 with r given in 

(ETJ). 

For the bosonic case, where we have (16T1) . this condition implies N + /(4N_) = 1 
and, accordingly, the stability is achieved by the ratio 

£ = 4. (81) 
Thus the redistributed numbers read 

N + = -N and AL = -N, (82) 
5 5 

which implies that a large portion (60%) of particles must be moved from one well to 

the other to achieve the equilibrium. As in the previous example, this portion does not 

depend on the particle number N. This second example also illustrates the fact that 

the bosonic case is macroscopically far from the mechanical equilibrium and hence the 

effect may easily be observed experimentally. 
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In contrast, the fermionic case has the limit (1661) . and the equilibrium condition 
r = 1 implies (N + /N_) 3 ■ (1 + 3/(2iV_)) _1 = 1. This is achieved with the ratio, 
N+ 1 1 

it " 1 + 2 F < 83) 
The outcome shows that the instability caused by the net force on the partition is 
extremely small so that even the redistribution by a single particle from W~ to W + can 
reverse the instability (AF(0) < 0). We thus learn that, in terms of redistribution, the 
effect of the quantum force on the partition is almost invisible for the fermionic case. 
Again, we see that the fermionic setting is macroscopically close to the mechanical 
equilibrium. 



7. Conclusion and Discussions 



In the present paper, we have investigated physical consequences of nontrivial boundary 
conditions in a quantum well realized by a partition placed at the centre of the 
well. We have examined the pressure appearing on the partition due to the distinct 
boundary conditions on the two sides (Dirichlet on the left and Neumann on the 
right), taking into account of the thermal effects under finite temperatures as well 
as the different statistics of the particles in the well. We have found that, for both 
bosons and fermions, the net force Af(t) acting on the partition in the quantum well 
is nonzero at the zero temperature limit t — > 0, and remains practically constant for 
extremely low temperatures before it starts to decrease gradually for higher (but still 
low) temperatures. Knowing that the energy spectrum is different in the two half 
wells W~ and W + separated by the partition for all energy regimes, this property 
is not unexpected. What is unexpected, however, is that this decrease stops at a 
certain temperature t m i n and afterwards the net force increases persistently up to the 
high temperature region t — > oo where it diverges according to the square root of the 
temperature y/i. 

We have also observed a salient scaling property in the particle number N for finite 
temperatures, that is, the force Af(t) is of the order of N for the bosonic case while 
it is of the order of N 2 for the fermionic case. Furthermore, the minimal force A/ min 
and the minimal temperature t m j n follow the same scaling law, i.e., they are both of 
the order of iV for bosons and iV 2 for fermions. For t — > oo, the difference in the 
statistics disappears and the force becomes proportional to N for both of the two cases. 
These properties are found numerically and have also been successfully reproduced by 
analytical approximations developed in this paper. 

This curious scaling dependence pertinent to the particle statistics may be 
understood intuitively by looking at the limit t — > 0. Namely, in the bosonic case 
the force A/(0) consists of the difference in the forces in the ground levels of the two 
half wells Af\ multiplied by the number N of the particles all of which reside in the 
same levels. In contrast, in the fermionic case the force A/(0) consists of the differences 
in the forces in the levels Af n = e~ — e+ = n 2 — (n — 1/2) 2 pa n (see ((SJ) and (|T3|) ) which 
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increase linearly up to the Fermi level n = N, thus yielding J2n=i n °^ in total. 
For finite temperatures, this observation will be valid until the statistical property of 
fermions becomes obscured for higher t where both of the forces tend to be proportional 
to N. For the quantum well in one dimension, this scaling dependence pertinent to the 
particle statistics may also be observed in other physical phenomena, not just in the 
quantum force we have studied. The shift of the partition mentioned in the previous 
section can be one such example. 

The difference in the scaling does not necessarily imply that the quantum force is 
more easily measured for the fermionic case than the bosonic case. In fact, we have seen 
in the (imaginative) shift in the partition or the redistribution of particles caused by 
the force that the effect of the force in these quantities is rather marked in the bosonic 
case while it is almost invisible for the fermionic case. 

Compared to the bosonic case, the fermionic case admits an additional property: 
for low temperatures the force exhibits a subtle step-like pattern, a small dent, just 
after it starts to decrease from the constant plateau in the extremely low temperature 
regime. This property has also been reproduced by analytical approximations, at least 
semiquantitatively. 

From the expectation that all the effects that arise from quantum boundary 
conditions should vanish at high temperatures where the classical picture would be 
available, the observed divergence of the force seems quite unusual. However, this 
may be understood by the fact that, contrary to most quantum systems in higher 
dimensions, one dimensional quantum wells possess energy spectra with increasing 
level spacing for higher energy levels (which is actually valid not only for wells with 
Dirichlet and/or Neumann boundary conditions but also for all other wells as well |15j). 
In other words, the dimensionality of quantum wells can be examined by their high- 
temperature behaviour, too. In actual realizations of the well, there is of course a 
maximal height of the potential, which will modify the high-temperature behaviour and 
eventually assimilate it to the classical one in the limit t — ► oo. 

It is also informative to consider what happens if the infinite well is replaced by a 
harmonic oscillator potential in one dimension. In the harmonic case, the levels follow 
each other with a constant spacing. This constant is the same on the two sides of the 
partition, but since the ground state energy is different the whole spectrum on one side 
is shifted by a constant with respect to the spectrum on the other. It turns out that the 
force difference for the nth levels on the two sides proves to be decreasing Af n « (7m) ~ 5 
for higher n as opposed to the linearly increasing behaviour Af n « n in the infinite well 
case. There are two sources for the different n-dependence in the harmonic case. One is 
that the energy level difference between the two half wells is constant and does not give 
larger contributions for higher n. The other is that because of the infinite stretch of the 
harmonic potential the higher energy states extend more in space towards infinity and 
are, correspondingly, less sensitive to a shift of the partition at the origin. As a result, 
one finds that the total net force decreases for high temperatures, with a t~^ asymptotics. 
In parallel, the iV-dependence also changes in the harmonic case and, for example at 
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low temperatures, it is of the order of A 7 "! for fermions (since J2n=i x while it 
remains A" for bosons. Therefore, the level spacing and the steepness of the confining 
potential both influence the various aspects of the net force considerably. The detail of 
the analysis for the harmonic potential case will be reported elsewhere. 

The calculation presented here could be performed for other, more general boundary 
conditions as well. For most boundary conditions, however, the energy levels will 
be determined by transcendental equations and, accordingly, additional technical 
difficulties, especially analytical ones, will arise. In this respect, the present combination 
by the Dirichlet and Neumann conditions is certainly special in that it admits a simple 
and yet distinct set of energy levels for the two half wells. Nevertheless, one may expect 
that, in this example where the two 'extreme' boundary conditions are used, most of the 
generic features that could arise from nontrivial boundary conditions have appeared. 

Our setting is also idealized from the point of view of possible direct experimental 
verifications. For instance, inspired by the second example mentioned in section 6, one 
may consider the measurement of the instantaneous current flow between the two half 
wells which is expected to occur when they are connected by a wire. Admittedly, at 
the present status of nanotechnology the boundary effects such as this may still be hard 
to be examined. However, the main message of this work is that the difference in the 
boundary condition creates mechanical and thermodynamical inequilibrium, and it is a 
realistic assumption that future technological developments will enable us to observe it 
in one form or another. We hope that the present study furnishes an intimate picture of 
the effect of boundary conditions in quantum mechanics, paving the way toward fuller 
understandings of quantum singularities in general. 
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Appendix A. Outline of the numerical calculation of the net force 

The purely numerical computation of the net force Af for a given particle number A" 
and temperature t used in this paper is performed in the following scheme. First, in each 
half well, we solve the total number condition (E]) for a, then substitute the obtained a 
into the corresponding expression (ITTj) for the force /, and at last we take the difference 
of the two forces /~ and / + . We have to control errors coming from three sources. 
The first is that we use finite truncations of the involved infinite sums, the second is 
that we obtain a via a numerical solution of the number condition, and the third is the 
floating-point errors of the calculations. 

For a prescribed precision of the force /, we estimate the required preciseness in a 
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by taking the partial derivative of the sum (1111) with respect to a at a fixed temperature 
parameter b. Similarly, the derivative of the sum (Q with respect to a is used to find the 
corresponding preciseness at which ([91) has to be fulfilled when seeking for the solution 
a. 

When setting this required precision, we also incorporate the error of the truncation 
of the infinite sums we face at, i.e., of Qj, ffTTT) . and their mentioned derivatives. To 
this end, simple upper estimate formulas for the dropped infinite terms are derived. For 
low temperatures, we can estimate from above by some appropriate geometric series 
q n or its generalization y]?L, (c 2 n 2 + cvn + c )q n , which is summable 

"trunc " — "trunc v ' ' 

in closed form. In parallel, for high temperatures, it is better to replace any dropped 
discrete sum with a corresponding continuous integral (cf. flB.ll) - flB.3l) ). Then we change 
the integrand to a simpler upper bound of it, to reach an integral like e~ y2 dy . 

Finally, we employ the standard approximate closed expressions to evaluate this simpler 
definite integral (the first few terms of its large-y tmnc expansion). 

The floating-point precision must also be chosen appropriately. The number of 
floating-point digits must be such high that a further increase in the number of digits 
influences the force difference result well within the uncertainty we prescribed for it. 
Unsurprisingly, we find that higher temperature t and larger particle number N require 
better floating-point precision. One should use such a mathematical computer software 
that allows to have 20 digits or more. Our calculations have been performed using Maple 
(© Maplesoft, Waterloo Maple Inc.). 

Appendix B. Calculations for the medium temperature regime 

For both the bosonic and the fermionic case, we will make use of an integral 
approximation of the infinite sums. Namely, the trapezoid approximation of integrals, 
applied for a function g(y) behaving 'peacefully' in [yi, oo) with lim^oo g(y) = 0, yields 



for constant intervals Ay = y n +i —y n = const. This approximation improves for smaller 
Ay. Remarkably, the trapezoid approach provides an approximation one order better 
(in Ay) than the simple rectangular one. 

Now, in all our applications used in the text, we have 




(B.l) 





(B.3) 



T-l/2 



in our case 
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In the main body of the present paper, wherever a sum is approximated by an integral, 
we have used (IB. 31) . 

Let us consider the bosonic case where we have t ~ N in the medium temperature 
region. There, we can show that the relative error of the approximation (126]) tends to 
zero in the large- N limit. Indeed, 



J_y, r^J 1 

N ^ L e a+ y™ - 1 a + 



^^^^^^^^^ 



' fl " 2r/ i-§] + ^/ I .... . -Kv 



N{ 4 V 6 J y/bJo \e a+ y -1 a + y 2 
]-\o(N ) +o(N 1/2 ) (B.4) 



iV 



x 



-1.2942 0.1842 

which is altogether an O (iV -1 / 2 ) quantity, vanishing in the limit N — > 00. 

Next, we derive ( 134)) and its order of relative error, under the conditions b ~ 1/iV <C 
1 and |a + |, |a~| ~ 1/iV <C 1 . For the force on one side, we can write 



°°p 00 « + ftp 00 z 

f — \ " hf 4- N — \ ^ - \ ^ n 

J ~ 2^ p a+be n _ 1 ' °i + iVa — 2^ p a+6e n _ 1 — 2^ £ n _ i 

n=l e 1 n=l e 1 n=l e 1 

with z n := a + be n . Thus, for the force difference, we have 

A(bf + No) = bAf + NAa = £ A (-^37) 



(B.5) 



n=l 



, , \ dz e z — 1 / 

n— 1 \ / z=z n 

Note that in the sum ( 1B.6I) . the factor 

Az n = z~ - z+ = Aa + bA(n - r) 2 = A{a + br 2 ) -2b Am (B.7) 

is small Az n = O (iV -1 / 2 ) up to n = O (i 1 ^ 2 ) , while the higher n terms in the sum are 
irrelevant as being exponentially suppressed. 

In ( ]B.6j) . the term f-j^ e *-i ) can er ther used for the quantity in W + or W~ , 
or even for the average of the two with increased preciseness. Fortunately, the two will 
prove to differ only in a subleading order, and we do not need to specify the choice. We 
proceed by rewriting Az n as 

Az n = [Aa + b(Ar 2 - 2tAt)] - 2ArVb y n (B.8) 

(a form that is valid on both sides W + and W~), with which we evaluate the sum ( IB. 61 ) 

as 



00 / d 7 

\Aa + b{Ar 2 - 2tAt)] V -f ~ 

' ' n=l " 

both terms ©(AT" 1 ) 



dz e 2 
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2ArVb J2 



z 



■ Vn 



V J I K J VbJ yi \dze z -l 



dy 



z=a+y* 



(B.9) 



0(6°) 



2ArVb 



O b°) + 



l r°° ( d z 



Vb 



in 



dz e 2 



■ydy 



z=ct-\-y 2 1 



dz 



Among the obtained 2+2 = 4 terms, the first three provide only an O \ N 1 ^ 2 ^j correction 
to the fourth term, which we calculate as 



At 



00 / d z 



dz e 2 



dz = At- 



e Zl - 1 
At [1 + O (a + yl) 

OiN- 1 ) 



Hence, we conclude that 



bAf + NAa w -- + O (N- 1/2 



or 



N 



1 t 

"2iV 



A(to) + C(iV- 1/2 ), 



+ C(iV~ 1 )(B.10) 

(B.H) 
(B.12) 



which is (1341). 
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